Preparation of stable excited states in an optical lattice via sudden quantum quench 
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We study how stable excited many-body states of the Bose-Hubbard model, including both the 
gas-like state for strongly attractive bosons and bound cluster state for repulsive bosons, can be 
produced with cold bosonic atoms in an one-dimensional optical lattice. Starting from the initial 
ground states of strongly interacting bosonic systems, we can achieve stable excited states of the 
systems with opposite interaction strength by suddenly switching the interaction to the opposite 
limit. By exactly solving dynamics of the Bose-Hubbard model, we demonstrate that the produced 
excited state can be a very stable dynamic state. This allows the experimental study of excited 
state properties of ultracold atoms system in optical lattices. 

PACS numbers: 37.10.Jk, 03.75.Lm, 03.75.Kk 



I. INTRODUCTION 

Recent experiments with ultracold atoms have of- 
fered exciting opportunities to study quantum many- 
body physics in a highly controlled manner [J-Q. Be- 
yond simulating the ground state (GS) properties of var- 
ious many-body systems, the uniqueness of cold atomic 
system, such as the low dissipation rate and the Fesh- 
bach resonance technique, has led to the experimental 
realization of stable excited states, as demonstrated by 
recent observation of a one-dimensional (ID) stable ex- 
cited state called super- Tonks-Girardeau (STG) gas @ 
and repulsively bound atom pairs [6]. In general, a sta- 
ble excited state is hard to be realized in traditional solid 
state systems since a pure excited state is not stable due 
to the energy dissipation between the system and the 
environment. The STG gas provides one of the coun- 
terintuitive examples realized in cold atom systems with 
no analog in solid state systems. It describes the lowest 
gas- like phase of the attractive Bose gas @, H| , which is 
however a highly excited state against its cluster-like GS. 
The stability of the STG gas could be understood from 
the dynamics of the ID integrable Bose gas [9]. An- 
other example, i.e., the repulsively bound atom pairs, is 
also counterintuitive at first glance since two atoms with 
strongly repulsive on-site interaction usually repel each 
other. Theoretical studies have revealed the existence of 
exotic repulsively bound pairs in optical lattices p^4l5| . 

The experimental realization of stable excited states 
provides a promising new area for searching novel quan- 
tum states in cold atom systems [161 Il7|. However, the 
diversity and mechanism of realizing stable excited many- 
body states are still not well understood. Questions arise 
whether the STG gas can also be realized in a lattice sys- 
tem as the Tonks-Girardeau (TG]_gas [1] and whether the 
two independent experiments [E|,l6| can be understood in 
the same theoretical framework. In this work, we study 



how to prepare specific stable excited states, including 
both the gas-like excited state in attractive regime and 
repulsively bound cluster of atoms, for a Bose gas in a 
ID optical lattice described by the basic Bose-Hubbard 
model (BHM). By exactly solving quench dynamics prob- 
lem of BHM, we first show that a stable gas-like phase 
for the Bose gas with strongly attractive on-site inter- 
action can be realized in the ID optical lattice by sud- 
denly switching interactions from the strongly repulsive 
regime to the attractive regime. Such an excited gas- 
like state avoids collapsing to the atom-cluster GS even 
under very strongly attractive interaction and could be 
viewed as a realization of STG gas in optical lattices. 
Furthermore, we show that a stable repulsively bound 
cluster state composed of iV-atoms can be also realized 
by sudden switch of the interactions from the attractive 
side to the repulsive side. The existence of the repul- 
sively bound cluster in optical lattice is closely related to 
the phenomenon of repulsively bound pairs. Our study 
suggests that the two seemingly unconnected phenom- 
ena p| |(| could be understood in a unified theoretical 
framework, i.e., they could be explained by calculating 
the sudden quench dynamics of the BHM with different 
initial states. 



II. SYSTEM AND SCHEME 

We consider the system of ultracold bosonic atoms in 
a ID deep optical lattice, which can be described by the 
BHM [11 III 



H 



•i+l 



1), (1) 



where b\ is the creation operator of bosons at the ith site, 
J and U denote the hopping strength and on-site interac- 
tion, respectively. The ratio U /J can be tuned by varying 
the depth of the optical lattice and using Feshbach reso- 
nance 18j. For convenience, we set J = 1 as the energy 
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scale. Despite the BHM having been studied by vari- 
ous methods [l9l - [2"lj ]. the model is generally not exactly 
solvable by analytical method [22[ . As the GS proper- 
ties have been extensively studied, the properties of ex- 
cited state and related non-equilibrium physics based on 
the Hubbard model have recently attracted lots of atten- 
tions [l|| [13, [23l - l27| . In this work, we propose a scheme 
of preparing stable highly excited states of the BHM via 
a sudden switch of interactions from the strongly repul- 
sive regime to the attractive regime and vice versa. Sim- 
ilar kind of quench has previously been used to create 
metastable states in other systems [H, H, [lj} ■ We shall 
demonstrate how an initially prepared GS translates to 
a highly excited state by solving the quantum dynamics 
of BHM. 

Suppose that the initial state \&i n i(t = 0)) = \ipo(Uo)) 
is prepared in the strongly interacting regime with either 
Uo > or Uq < 0, after a sudden switch to the opposite 
regime with interaction strength U , the wave-function 
= e -iH{U)t \^ mi {jj )) can be calculated via 
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where c„ = {ip n {U)\ ipo(Uo)} with \ip n (U)) representing 
the n-th eigenstate of the BHM with on-site interaction 



U. It is straightforward that 



is the transition prob- 



ability from the initial state to the n-th eigenstate of 
H(U). To study the quench dynamics of the BHM, we 
shall scrutinize the full spectra and eigenstates of the 
BHM both analytically and numerically. 
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FIG. 1: (Color online) Full spectra for the two-particle BHM 
with L = 50, U = 10 (left) and U = -10 (right) in K- 
space. Left inset: The transition probability from the initial 
scattering GS with Uo = 80 to the lowest scattering state 
(P 3 s) and the dimer GS (P s b) for various attractive U. Right 
inset: The transition probability from the initial attractively 
dimer GS with Uo = —80 to the repulsively dimer state {Pbb) 
and the scattering GS (Pbs) for various repulsive U. 



The complete solutions of the B AE are composed of C£ 
real solutions (scattering states) and L complex solutions 
(bound states). The former ones follow by solving the 
logarithm of the above BAE: 



III. BETHE-ANSATZ SOLUTION FOR 
TWO-PARTICLE BHM 
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We begin with the two-particle problem of the BHM 
which is exactly solvable with the aid of Bethe-ansatz 
(BA) method. Although the two-particle problem is 
quite simple, its analytical result can provide us quite in- 
structive understanding to many-particle systems as we 
shall discuss later. The BA wavefunction takes the form 
of I*} =Y, XuX2 ^{xi^2)\xi,x 2 ) with 



^(x 1 ,x 2 ) = A 12 e 



i(kixi+k 2 x 2 ) 



+ A 21 e 



i(k 2 xi+kix 2 ) 
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defined in the domain x\ < x 2 . The wavefunction in the 
other region x 2 < x\ can be obtained by the symmetry 
of wavefunction. Explicitly, the coefficients are given by 
Ap = (— l)- p (sinfc P2 — sin/c pi —ill/2), where P = {pi,p 2 } 
is one of the permutations of 1, 2 and (— l) p = 1 or — 1 for 
even or odd permutation. Under periodic boundary con- 
dition, the quasimomenta kj (j = 1,2) fulfill the Bethe- 
ansatz equations (BAE) 



exp (ikjL) = 



sin ki — sin kj — iU /2 
sin ki — sin kj + iU /2 



The total momenta of the system is given by K = k\ + k 2 



and the eigenenergies are given by E = — 2^iLi cos/ci. 



with/j = -(i-l)/2, ...,(£- 1)/2 and h < h- While the 
bound states correspond to the string solution with the 
form of k\ = k + iA and k 2 — k — iA, where k and A are 
real. In this case the total momenta and eigenenergy take 
the form of K = 2k and E — — 4 J cos k cosh A, respec- 
tively. In terms of k and A the original BAE transform 
into finding root of 



cos (kL) exp (AL) 



2cos/csinh A - U/2 
2 cos £; sum A + [7/2' 
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where k = Ijtt/L with Ij = 0, ±1, ±2, ±L/2 - 1. Be- 
sides these L — 1 bound states the Lth one corresponds 
to Kl — n and El = U/J. We obtain all eigenstates 
by solving (j4]) and ([5]). Full spectra for example systems 
with U = ±10 and L = 50 are displayed in Fig. [TJ 

Taking advantage of analytical results for the two- 
particle case, we can calculate the transition probabilities 
exactly. First, we consider the case in which the initial 
state is prepared as the GS of the BHM in the repulsive 
side. After suddenly switching the on-site interaction 
to the attractive regime, we calculate the overlap of the 
initial wavefunction with the eigenstates of the attractive 
BHM. In the left inset, we show the transition probability 
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from an initial GS (marked by the cross in the left Fig.l) 
with Uq — 80, to the lowest scattering state (marked by 
the cross in the right Fig.l) with different values U < . 
Here the transition probability from the initial scattering 
GS to the final lowest scattering state is denoted by P ss 
which is very close to one after switching to the strongly 
attractive regime. On the other hand, the probability 
for dynamically falling into the attractively bound state 
(denoted by P s b) is almost zero. In the right inset, it 
is shown that after switching to the strongly repulsive 
regime, the initial attractively pair state (marked by the 
star in the right Fig.l) is transformed to the repulsively 
bound state (marked by the star in the left Fig.l) with 
the transition probability close to 1. Here the symbols 
Pbs and Pbb represent the transition probabilities from 
the initial bound state to final scattering state and the 
initial bound state to final bound state, respectively. 
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IV. MANY-PARTICLE SYSTEM 

For the many-particle system with N > 3, the BHM is 
no longer exactly solvable by the BA method which can 
not properly treat the multi-occupation case [22| . Never- 
theless, for a finite size system, we can resort to the full 
exact diagonalization (ED) method to calculate the full 
energy spectra and eigenstates. Consequently the tran- 
sition probabilities from the initial GS to arbitrary final 
states are straightforward to be calculated. In general, 
it is a formidable task to get the full spectra of a large 
system as the basis dimension of a iV-particle BHM with 
size L is given by D = (N + L - 1)1/ [N\ (L - 1)!]. 

Despite the full spectra becoming very complicated as 
the particle number increases, we can still find some com- 
mon characteristics of the spectra for systems with differ- 
ent sizes. When \U\ 3> J, the spectra is split into a series 
of separated bands. For the repulsive case, the lowest 
band is a scattering continuum of N asymptotically free 
particles, whereas the top band is a narrow band formed 
by the iV-particle repulsively bound state. Between the 
top and bottom bands, there exists a series of scattering 
continuums formed by bound cluster states or formed by 
bound state and free particle. To give a concrete ex- 
ample which may guide us to understand the structure 
of the spectra, we display the energy-momentum spectra 
for a system with N — A, L — 30 and U = ±15. As 
shown in Figj2^,, the bound states on the top band have 
the energies about QU , and the lowest band is a scat- 
tering continuum of 4 free particles. In between, the en- 
ergies of three separated bands are approximately given 
by 3U, 2U and U, which suggests that these bands cor- 
respond to, respectively, the scattering continuum of a 
trimcr and a free particle, two dimers, and a dimer and 
two free particles. The spectra for system with attractive 
interactions has a similar structure in reverse order. The 
zero-momentum iV-particle attractively bound state is 
the GS, whereas the scattering continuum of N free par- 
ticles is on the top of the spectra corresponding to highly 



FIG. 2: (Color online) Full spectra of energies vs momentum 
for the BHM with N = 4, L = 30 and U = ±15. The symbols 
of cross represent the lowest scattering state in the bottom 
band (GS) for the repulsive U (left) and the lowest scattering 
state in the top band for the attractive U (right), whereas 
the symbols of star represent the highest repulsively bound 
cluster state (left) and the lowest attractively bound cluster 
state (right). 

excited states of the attractive system. We note that 
these separated bands are no longer discernable when the 
interaction strength is comparable to the band width. 

For a sudden switch of the interaction from the 
strongly repulsive regime to the attractive side, we evalu- 
ate the transition probabilities from the initial repulsively 
GS to the final states. As displayed in Figj3k, the tran- 
sition probability to the lowest excited state in the top 
scattering band is very close to 1 in the strongly attrac- 
tive regime, whereas the probability for falling into the 
attractively cluster GS is almost completely suppressed. 
Since the transition rate to the scattering phase is very 
close to 1 in the strongly interacting regime, we expect 
that one can prepare such a highly excited state exper- 
imentally through switching the interaction from strong 
repulsion into strong attraction following the same way in 
the experiment of the STG gas. Actually, the stable ex- 
cited scattering state prepared in this way can be viewed 
as a realization of STG gas in optical lattices. The gas- 
like excited state is no longer stable and decays quickly 
if the system enters to the weakly interacting regime. 

We can understand the stability of the lowest scatter- 
ing state in the strongly attractive limit from the ana- 
lytical BA solution. We note that the BAE solutions 
for the repulsive case and that for the STG gas corre- 
spond to the same set of {/,■} according to d3J). In the 
limit \U\ — > oo, the solutions given by kj — Ij2n/L 
are exactly the same, which means that the repulsive 
GS state and the lowest scattering state in the attrac- 
tive side are identical in the infinitely interacting limit. 
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FIG. 3: (Color Online) (a) The transition probability from 
the initial GS of the repulsive system with Uo = 100 to the 
lowest gas-like state (Pas) and the cluster GS (P 3 b) after the 
switch to attractive side, (b) The transition probability from 
the initial GS of the same system with Uo = —100 to the 
GS (Pbs) and the repulsively cluster state (Pbb) for various 
repulsive U. 



In the strongly interacting regime with \U\ 3> 1, the 
quasi-momentum distributions for the repulsive TG gas 
and the STG gas approach the free fermion distribution 
from different sides, and consequently the overlap be- 
tween the repulsive GS and the STG state is close to 
1, i.e., \(^stg(-U)\MU))\ -»• 1 as \U\ -> oo. The 
above exact discussion for 2-particle system can be di- 
rectly extended to the many-particle systems. This is 
based on the observation that, in the strongly interact- 
ing regime, the extended BA solutions ^ can be used to 
describe the properties of scattering states very precisely 
for both the repulsive and attractive cases, although 
the solution is not an exactly analytical solution for the 
many-particle system in the rigorously integrable mean- 
ing |22j. Here the many-body BA wavefunction takes 



the form of ^ (x\ 
coefficients Ap - 



, acjv) = Yip Ape 1 ^? kp j Xj , where the 



P = {pi,p2, ... 7 pn} is one of the permutations of 1, N, 
and J2p is the sum of all permutations. The quasimo- 
menta are determined by the BAE 



kjL = 2nlj 
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Both the GS solution for U > and the STG solution for 
U < are determined by the same set of {Ij} = { — (N — 
l)/2, (N + l)/2}. To confirm that, we calculate and 
compare the GS energy for the U > case and energy of 
the lowest scattering state for the U < case by both the 
BA solution and ED (see FigfJJi). The energies obtained 
by the two methods agree very well with, for example, 
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FIG. 4: (Color Online) (a) Etg and Estg vs \U\ for sys- 
tem with N = 6 and L — 10. Dash and solid lines are 
obtained by solving BAE whereas the squares (STG) and 
dots (TG) are obtained by ED. (b) 6 vs \U\ with S(U) = 
\E E d(U) - Eba(U)\/\E E d(U)\. Comparison between the B A 
solutions and ED results suggests that the BA solutions can 
give very good results consist with the ED results in the 
strongly interacting regime, (c) The local two-particle cor- 
relation function Q2(U) — dE(U)/dU vs interaction. 



<S = \E E d - E B a\/\Eed\ < 10~ 5 for \U\ = 100 as shown 
in the FigHJs. As \U\ — > oo, we have Etg = Estg- 
Such a gas-like state shares similar characteristics with its 
continuum correspondence. For example, the STG phase 
exhibits stronger local correlation than the repulsive Bosc 
gas as displayed in the FigHb. 

Finally, we turn to the case in which the initial state is 
chosen to be the GS of the BHM in the strongly attrac- 
tive interaction limit, i.e., the cluster- type bound state 
with all the atoms tending to stay together due to the 
strongly on-site attractive interaction. After the on-site 
interaction is suddenly switched to the strongly repulsive 
regime, we also find that the cluster-type bound state is 
stable and the initial system translates to the repulsively 
cluster state as marked by star in Figf5^. In FigJ5]D, we 
show the transition probabilities from the initial attrac- 
tively bound state to the repulsively bound state and the 
scattering GS in the repulsive side. We find that the tran- 
sition probability to the repulsively bound state is close 
to 1 whereas the transition to the scattering GS is almost 
completely suppressed in a wide range of interaction. 
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V. SUMMARY 

In summary, we have studied the transition from the 
GS of strongly repulsive or attractive bosons to the gas- 
like highly excited state of attractive bosons or the repul- 
sively bound state through a switch of interaction. By 
calculating the transition probabilities, we have shown 
that the gas-like excited state and repulsively bound state 
are stable in the strongly interacting regime and thus are 
possible to be observed with cold atoms in optical lat- 
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